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I. INTRODUCTION 

Rapid progress has been made recently in computations of the free energy density, at 
nonzero temperature (T), of massless g'^cj)^ theory |jl[|, quantum electrodynamics (QED) 
and quantum chromodynamics (QCD) [Q, to three-loop (fourth order in coupling). 
For QED, the fifth (e^) order contribution has also been obtained 0] by dressing the 
photon lines of the three-loop diagrams. 

Compared to QED, a fifth order calculation in QCD will be more involved because 
gluonic self-interactions imply that many more lines in any three-loop diagram can be 
soft (i.e. at zero Matsubara frequency in the imaginary time formalism), and so must be 
dressed in order to obtain the full g^ contribution. In this respect a fifth order calculation 
in QCD will resemble the same order calculation in g'^ip'^ theory. In this paper we will 
compute the order g^ T^ contribution to the pressure of massless g'^cj)^ theory and show 
that by an optimal choice of rearranged Lagrangian it is possible to obtain analytical 
results without too much effort. 

The theory we are concerned with is defined by the Euclidean Lagrangian 

C, = \{d<pf + '^c^\ (1) 

where /z is the mass scale of dimensional regularisation((i = 4 — 2e) . We use the 
imaginary time formalism in which the energies take discrete values, 27rnT, n E Z. For 
perturbative calculations beyond leading order, it is necessary to take into account in 
a systematic manner the non- negligible collective effects 0. For the theory defined by 
eq.(|l]) this means that one uses instead the shifted Lagrangian ^j 

C, = [Co + ImV) - \m^' (2) 

with IV? = g'^ T^/24 the thermal mass generated at one-loop. The term within brackets 
in eq. (^ defines a dressed propagator while the last term is a new two-point vertex which 



prevents overcounting. With eq. (0) one can proceed to calculate any Green's function in 
the theory, order by order, in a consistent way. However for the calculation of Green's 
functions with static (zero energy) external legs, eq.(|D is not very economical since it 
involves some extraneous resummation. 

Recall (see, for example, 0,H) that for a static Green's function, its physical definition 
is already given in imaginary time, without the need to analytically continue to real 
time. Thus the power counting of infrared (IR) divergences may be safely done using 
the Euclidean propagators with discrete energies. Then, only the propagators at zero 
Matsubara frequency do not have an IR cut-off of order T, and it is only for these zero 
modes that the thermal mass ~ (7 T is a relevant infrared cut-off. Thus instead of eq. (0) 
one can use 

C= (Co + -mV^f^po.o j - ■^m^4>'^Spo,o (3) 

for static calculations, with po the energy in Fourier space of (f){x). The Lagrangian in 
eq.(H) is precisely what is suggested by the Braaten-Pisarski resummation scheme 
whiu involves dressing 'soft' lines with 'hard-thermal loops'. For the calculation of static 
quantities in imaginary time, the only soft-line is the zero-mode propagator and the only 
hard-thermal loop is the static one-loop thermal mass. 

We remark that a resummation as in eq.(^), which involves dressing only the zero- 
modes, was used very effectively by Arnold and co-workers 0,^ for their free energy 
calculations in gauge theories, but they used the more general expression (0) (as in Ref. 
|l[]) for their scalar free-energy calculation. The advantage of using minimal resummation 
(0) also for the scalar case can be seen by the following example. The propagator from 
eq.(l) is 

A(K) ^ iijliMl + ^^0-0 . (4) 
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Then the one-loop integral 
where we have used the notation 

Ko(even) /co(even) ^ ' 

with K^ = {ko,k). Notice that the integrals in eq.(^) were computable in closed form 
and the result is a term of order g^ and another of order g. Since the bosonic one-loop 
integral @ and others related to it occur frequently in Figs.l(a-g) , one easily sees that 
all those diagrams may be evaluated easily. By contrast the propagator of eq.(^ is 
1/{K'^ + m'^) and its one-loop integral is only available as a high temperature expansion 
0, making the evaluation of the diagrams of theory (0) more involved. 

Another point that needs clarifying is the choice of m^ in eq.@. In four dimensions, 
rn^ = ^^2A~- However since we are using dimensional regularisation, one might wonder 
if the value of m? in ^-dimensions should be used. For example, in |]l| the value m? = 
^-^2A — 1^"^^ ^^^ used (it being the hard thermal loop in 4 — 2e dimensions), while in 
0] the full one-loop value in (i-dimensions was used. Will these different choices affect 
the final result? If one is interested in renormalised values as e — > 0, the answer is 
no. Here is the proof: Let the resummation in eq.(|^) be done using some m'^{e) with 
'^^(O) ~ ^'2^- Now keep the e - dependence of m^{e) implicit, even when it hits 1/e terms. 
Since the full Lagrangian in eq.(0) is massless, no ultra-violet mass renormalisation is 
needed in perturbative calculations using dimensional regularisation. Therefore any "^ ^'^' 
terms generated when calculating a renormalised quantity (i.e. after including coupling- 
constant renormalisation) must mutually cancel. Finally only terms of the form m?{e)e"' 
{n > 0) will appear and then, as e — >■ 0, only m'^{0) survives. End of proof. 



In explicit calculations one finds that if the e-dependence of m^(e) is expanded 
out then extra finite (e'') terms are generated from some diagrams because of the e- 
dependence of the mass in the propagator, but these will cancel , order by order in g, 
with similar terms generated from e-dependence of the thermal counterterm. The fact 
that one can use m^(0) rather than some m^(e) clearly simplifies calculations and avoids 
redundant cancellations. On the other hand one can exploit the proof in the last para- 
graph to provide a cross-check on calculations. That is, one can use m^ = ^^^(1 + e^), 
where A is an arbitrary regular function of e ( which may also depend on g and T/fi), 
in eq.(|^) and verify that one's final renormalised result is independent of A as e ^ 0. 
Such a check had been performed P] for the pole of the propagator of g'^(j)'^ theory to 
two-loop 0, and we have also done it for the calculations in this paper. However in 
order to emphasize the simplicity of using ?7i^(0), we will present here the results for this 
case only. 

In the next section we reconsider the calculations of the pressure to order g^ using 
eq.(^) and recover the results of |jl],0. Then in section-Ill the g^ terms are obtained and 
the results are summarised in section- IV. We conclude with some comments in section-V. 

II. LOWER ORDERS 

The diagrams which contribute to the pressure to order g^ are shown in Fig.(l). The 
g'^ terms from these diagrams were extracted in refs. |jl|,^. Here we reconsider those g^ 
terms using the minimal Lagrangian, eq.(|), with m^ = ^^-- The reader comparing the 
results here with those of 0,^ should note this difference; order by order (rather than 
diagram by diagram) our results here agree with those of [|I|,Q . Diagram(a) contributes 



Va = -\l[dP]\nA-\P) 



-i JldP] ln{p2(l - <5,„,o) + iP' + m')5,,,o} 



(8) 






^2X(^-l)/2 
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In eqs.(§,P), the first term represents the ideal gas contribution while second term is the 
plasmon contribution (in dimensional regularisation, see [0) obtained by dressing the 
zero-mode of the one-loop diagram. 

The contributions of diagrams (b) through (g) are: 
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Note that in Vd only the one-loop ultra-violet (UV) coupling constant renormalisation 
counterterm is used. (We are also using minimal subtraction) . The g'^/e piece is required 
to cancel similar pieces from diagram (g) and (h). The g^/e and g^/e divergences in 
eq.(|lTD are due to the mixing of IR resummation effect with the UV renormalisation 
and will cancel against similar terms generated from diagrams (g) and (h). 

The order g"^ contribution from diagram (h) is obtained by setting m = in the 
propagators since the integrals are IR finite. The result Vha has been evaluated analyt- 
ically by Arnold and Zhai and we simply quote their value (which agrees with earlier 
semi-analytical evaluations in |1|,0|), 

^ ^^^-^^^I^l,7 + 181n^-^-12^p^ + 48^p^-187 + — j . (13) 

Thus the sum of diagrams up to order g"^ is 
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The terms to fourth order agree with 



III. FIFTH ORDER 

We now pick up the subleading pieces from diagram (h). For this we first rewrite 
A{K), given in eq.(|), as (c.f. Ref. i) 



A(P) 



rn? 5, 



Po,0 
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= Ao(P) + A*(P). 



(15) 



Then 



^ '" 48 
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where 



lo = //u'^ JldKdQ dP]Ao{K)Ao{Q)Ao{P)Ao{K + Q + P) 
Ii=g^fi''J[dKdQdP]A*{K)Ao{Q)Ao{P)Ao{K + Q + P) 
h = (7V' j[dKdQ dP]A*{K)A*{Q)Ao{P)Ao{K + Q + P) 
h = /Ai'' j[dKdQ dP]A*{K)A*{Q)A*{P)AoiK + Q + P) 
I^ = g^fi^' f[dKdQdP]A*{K)A*{Q)A*{P)A*{K + Q + P) 

The integral Iq contributes to Vh^ and was considered in the last section. We will now 
extract the order g^, g^/e and g^/e pieces from Ji through I4. Though our final objective 
is to calculate the pressure only to fifth order, we have to ensure that all subleading 
divergences such as g^/e (see eq.([Tl])) cancel (there are no divergences beyond g^/e from 
these diagrams). 
Consider 
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Since the external coefficient is 0{g^), we need only the order g piece from the {Q, P) 
integrals. Write the (Q, P) integrals as (the following discussion parallels that of I sun in 
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As the second term above is IR safe, one can expand the denominator 
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so that the second integral in eq.(|T9D is 

MQ^^It^^Jtt^ -2m l[dQdP]^-^f^^^l.{q + p) + 0{m') . (21) 
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The first term in eq. ( pi]) vanishes in DR 0,@,|I while the second term vanishes when the 
final k integrals are performed in eq.(|l8]). Hence 

[dk) f {dpdq) 



^--^vw-^V^,|^/. 
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F(/c2 + l)7 p-iq^lq + p + mky 
The {p, q) integrals are logarithmically sensitive to m in the infrared. They also have a 
logarithmic UV singularity in d—1 = 3 dimensions and so must be evaluated in c? = 4— 2e 
dimensions. After a standard evaluation of the {q,p) integrals one can perform the final 
k integrals easily by keeping only the terms to 0(e°). We obtain 
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Next consider I2: 

T =.4 4 6.7.2 f jdkdq) r [dP]5ko,oSg,^o 

' ^ ^ J k^q^{k^ + m^){q^ + m^)J P\K + Q + P)^ 
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The P-integral is 

J p^[j)^ mk + mqY J -< pg + (p + ?7i/e + ?Tigj "^ 

Thus 

+ /.^^,^m-^^T3 / ^^^H^L^ ^ 0(7) , (24) 

The first hne in (|^) is of order 0{g^). Since we require at most the 0{g^/e) piece, we 
can set c? = 4 everywhere there except in the P-integral which gives the pole 

m'-^ = J^^+0(.°), (25) 



and so the first line of eq.(p4D is 



„4 2Tn2 

|^ + 0(A°). (26) 
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The integral in the second line of eq.(24) is finite as e — > 0. Therefore we need to evaluate 

1 f d^kd^pd^q 



(27r)9y k\k^ + l)p\p^ + l)q^{k + p + q)^ 
We decouple the {k,p, q) integrals by writing 

1 _ /"^3 ^^{k+p + q + w) 



(k + p + q)^ -I w^ 



m;2 J (27r)3 
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Inserting eq.(P^) into eg. (P7|) , the {w,k,'p,q) integrals become trivial, giving 

eq.(27) 



1 
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We have defined a function 



°° dr 
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which will appear repeatedly. The use of eq. (^8]) in (^) is equivalent to evaluating the 

(see 0). In ( pOf ) one recognises ^ 



47rr 



momentum integrals ( |27D in coordinate space ( 

as the coordinate space Coulomb propagator and |^ as the coordinate space screened- 

Coulomb propagator. Adding eqs. (p6D and (0) gives 
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Finally, Is and I4 are both finite and their evaluation is analogous to the steps leading 
from eq.(^) to ( ^Of ) and utihses eq. (|3TD . We find 
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The full g^ contribution is then obtained by adding this to the value of diagrams (a-g), 
some of which also contain g^ pieces, given in eqns.(^|l2|). The result is displayed in the 
following section. 



IV. SUMMARY OF RESULTS 
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As required, all divergences, including the spurious g^ /e terms, have cancelled. The 
renormalisation scale /i appears explicitly in the In (^ j terms and also implicitly in the 
coupling constant, g = g{fj.) . One can eliminate the In ( ^j terms by re-expressing the 
pressure in terms of the temperature dependent coupling g(T) given by 
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A simpler but perhaps less instructive way to obtain eq.(^) is to choose /i = T in 
eq.®. 

V. CONCLUSION 

Using the minimally rearranged Lagrangian (§) together with an e-independent ther- 
mal mass m^ = ^^^, we have verified previous |l|J^ fourth order results for the pressure 
of massless g'^(j)'^ theory and then extended the calculations to fifth order. Our final 
result is given by eqs. (|36D and (^). At the fifth order a coupling constant logarithm 
appears for the first time in the pressure oi g'^(f)'^ theory. For QCD such coupling constant 
logarithms appear already at fourth order, but they do not occur in the pressure of QED 
(at zero chemical potential) because there are no self-interactions of photons (the only 
soft fields in imaginary time) and the conclusion follows by power counting |Q. 

It is natural to contemplate next a fifth order calculation in QCD. Based on our 
experience with g^(l)^ theory and QED we expect such a calculation to be technically 
simpler than the corresponding fourth order (three-loop) calculation: Nontrivial three- 
loop diagrams which are IR finite in the bare theory and computationally difficult (e.g. /q 
in g'^cj)^ theory), contribute at subleading order (g^) when at least one of bare propagators 
is replaced by a zero- mode dressed propagator (c.f. Ji in eq.(p!7D) so that the sum- 
integral over the dressed momentum line collapses to a three-dimensional UV finite 
integral, leaving only two overlapping frequency sums at most. Since frequency sums 
are the main complication in these calculations, this reduction saves effort. In practice 
further simplification has been observed : for the fifth order QED calculation, summing 
over gauge-invariant sets of diagrams results in the cancellation |Q of the terms with 
two overlapping frequency sums; for the scalar calculation in this paper the terms with 
two-overlapping frequency sums were found to contribute at higher order (g^) (see the 
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evaluation of Ji). That is, both the QED and scalar fifth order calculations turned out 
to be easier than expected. This bonus might prevail for QCD. 

As noted by Linde many years ago, the perturbative evaluation of the pressure 
in QCD breaks down at order g^ because of the absence of magnetic screening at lowest 



order. Braaten [10| has recently proposed a solution whereby one can obtain the coeffi- 
cient of the g^ contribution as a functional integral in a dimensionally reduced effective 
theory obtained by integrating out the hard fields in the original QCD path integral. 
Braaten has also described how his effective Lagrangian may be used to obtain the lower 
order g^ term and it would be interesting to compare the result of that approach with 
one using a shifted Lagrangian a la eq.(3). 

Figure Caption 

Fig.l(a-h) : Diagrams which contribute to fifth order. The propagators are given by 
eq.(4), the cross represents the thermal counterterm of eq.(3) while the blob in dia- 
gram(d) is the ultraviolet counterterm. 
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